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The prospects for high thermoelectric performance in phosphide skutterudites are investigated based 
on first principles calculations. We find that stoichiometric C0P3 differs from the corresponding 
arsenide and antimonide in that it is metallic. As such the band structure must be modified if high 
thermopowers are to be achieved. In analogy to the antimonides it is expected that this may be 
done by filling with La. Calculations for LaFe^Pn show that a gap can in fact be opened by La 
filling, but that the valence band is too light to yield reasonable p-type thermopowers at appropriate 
carrier densities; n-type La filled material may be more favorable. 



-a 
c 

o 
o 



> 

m 
O 

-1— > 



1 

C 

o 
o 



X 



There has been considerable recent interest in the elec- 
tronic and thermal transport properties of skutterudites. 
This is driven primarily by the discovery of two new 
high ngrformance thermoelectric (TE) materials in this 
class.Era The TE performance is characterized by a di- 
mensionless figure of merit ZT — aS 2 T/k, where a is 
the electrical conductivity, S is the thermopower and k 
is the thermal conductivity; ZT up to 1.4 at T = 600 K 
has been measured in skutterudites. Much of the effort 
has focused on antimonidesB based on the expectation of 
lower values of the lattice thermal conductivity related 
to the heavier mass atoms as well as the likelihood of 
better carrier mobilities due to the chemistry of Sb as 
compared to say P. In fact, the two hifth ZT compo- 
sitions discovered are both antimonidesalj: CeFe^Sb^ 
and La{Fe, Co)iSbi2. 

The high ZT values in these compounds derive from 
two important features: (1) high power factors aS 2 re- 
lated to their particular electronic structures, which are 
apparently different both between the two compounds 
and from the corresponding binary, CoSb^; and (2) a 
strong suppression of the thermal conductivity of the bi- 
nary upon filling. This latter effect, though crucial for 
the TE performance, is understood only qualitatively in 
terms of phonon scattering related to rare earth vibra- 
tions. Attempts to obtain even better performance by 
various alloying and substitutions on each of the three 
sites have thus far been unsuccessful, although there are 
still many possibilities remaining to be explored. These 
efforts are complicated by the large variety of realizable 
modifications of these skutterudites and the general lack 
of detailed understanding of their effects on properties 
relevant to TE. Moreover, TE performance typically is a 
strong function of the doping level, further complicating 
the search. 

In the present brief report we present electronic struc- 
ture calculations for the phosphides, C0P3 and LaF e^Pyi 
and discuss these in terms of the implications for TE 
performance and in relation to the corresponding anti- 
monide materials in order to elucidate trends. All previ- 



ous first principles calculations point to a particularly im- 
portant role for bands associated with the chemical bond- 
ing of the pnictogen 4-membered rings in the skutterudite 
structure in determining transport properties. - a point 
that was emphasized early on by Jung et al.Q based on 
tight binding calculations. .Calculations for CeFe^P^, 
C eF e\Asn and CeFe^Sbi^u have shown these materi- 
als to be hybridization gap semiconductors with decreas- 
ing gaps as the lattice parameter increases and Ce-f hy- 
bridization decreases down the pnictogen column. Previ- 
ous calculations for the binaries C oSb^ And C oAs^ reveal 
generally similar electronic structures Jj but with differ- 
ences that are particularly significant in the region near 
the Fermi energy (E^) that dominates electronic trans- 
port. CoSbs is a narrow gap semiconductor with a highly 
non-parabolic valence band dispersion, while CoAs^ was 
found to be a zero gap semiconductor with parabolic 
bands. 

Zhukov has reported first principles band structure cal- 
culations for C0P3 finding-toe material to be a narrow 
indirect gap semiconductorulij. The relatively heavy con- 
duction bands with their multi- valley minima would seem 
initially favorable for the electronic aspect of TE per- 
formance with n-type doping. However, the calculations 
were done with the linear muffin tin orbital atomic sphere 
approximation method (LMTO-ASA). Because the skut- 
terudite crystal structure features large voids, low site 
symmetries and strong covalent bonding, such calcula- 
tions are particularly difficult, and in such cases may have 
band shifts of several tenths of an eV compared to more 
accurate general potential calculations. Because of the 
small indirect gap, this is enough to qualitatively change 
the picture from a transport point of view, implying the 
need for a general potential investigation as presented 
here. 

Our calculations were done in the framework of density 
functional theory using the general potential linearized 
augmented plane wave (LAPW) methodlJ which docs 
not make any shape approximations and uses a flexi- 
ble basis set including LAPW functions and local or- 
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bital extensions^ to relax linearization errors and treat 
semicore-states. Valence states were done in a scalar rel- 
ativistic scheme while fully relativistic calculations were 
done for core states in the atomic spheres (Rmt =2.1, 
2.1, 1.9, 2.5 a.u. for Co, Fe, P and La respectively). The 
the basis set convergence was tested using R m in ■ ^max 
from 5.0 to 8.5; a value of 7.0 was found to yield a rea- 
sonable computational effort with only a small error with 
respect the highest R mm ■ k MAX (AE T = 0.3 Sgjg). We 
used a (4,4,4) special points grid for the Brillouin zone 
integration, which we found to be converged. The elec- 
tronic density of states (DOS) was based on a 35 k points 
tctrahcdral sampling in the irreducible BZ . The Hcdin- 
Lundqvist parameterization for the exchange-correlation 
LDA functional is used. 

As discussed below, we find, for C0P3, a globally simi- 
lar band structure to Zhukov, but with changes near Ep 
that are large enough to drastically change the picture in 
an unfavorable direction from the point of view of TE. 

As mentioned, the binary CoSb 3 has a relatively high 
thermal conductivity and a highly non-parabolic valence 
band dispersion, which is unfavorable for high p-type TE 
performance. Meanwhile La(Fe, Co)^Sb\2 has both a 
strongly reduced k and a band structure that is modi- 
fied in such a way as to improve the electronic properties 
by shifting the valence band edge downwards due to Ex- 
pulsion from the La /-resonance above the Fermi levelll3. 
One may conjecture that a similar effect could be present 
in La(Fe, Co)^Pi2 as the band edge states have the same 
character, and if so the question arises as to whether the 
electronic properties relevant to TE may be improved. 
There is also interest in the electronic structure of filled 
phosphide skutterudites because of the observation of su- 
perconductivity in LaFe^P\2 with critical temperature, 
T c = 4.1 KB. 

The skutterudite structure (space group im3) consists 
of a simple cubic transition metal sub-lattice partially 
filled by almost square pnictogen groups (P4). Three 
quarters of these sites are filled with such rings oriented 
in [100], [010] and [001] directions according to the cubic 
symmetry and the remaining one quarter are left empty. 
In filled skutterudites these remaining sites are occupied 
by a rare earth ion, which modifies thermal and elec- 
tronic properties. Two symmetry independent parame- 
ters u and v determine the position of the P with respect 
to the metal ion in the center of the cubic cell; they con- 
trol the size and the squareness of the rings. We start 
from experimental crystal structure fixing the position 
of pnictonen group with respect the transition metal to 
u e = 0.1453 a and v e — 0.3482 a where the lattice param- 
eter is a = 0.77073 nm for C0P3 and u e — 0.1504 a_and 
v e = 0.3539 a where a = 0.78316 nm for LaFe^Pi^B. 

The band structure of C0P3 , as given in Fig.|l|, is metal- 
lic due to the fact that the pseudogap near Ep, which is 
characterstic of skutterudites, is entirely crossed by a sin- 
gle mostly phosphorus p band. This is the same band as 
crosses the pseudogap in CoSb 3 and CoAs 3 , but in C0P3 
it crosses the conduction bands above the Fermi level. As 



such, C0P3 is not very promising for TE applications un- 
less filling or other modifications alter the band structure 
enough to open a gap. The corresponding DOS and pro- 
jections are shown in Fig.||. In the most relevant region 
for TE properties (near Ep) our LAPW band structure 
is quite different from the LMTO-ASA resultsE3 which 
predict an indirect T-H energy gap. As already found in 
Ref. |§| for antimonides, there are (1) a single degenerate 
band, (2) a two-fold degenerate band and (3) a three- 
fold degenerate band at T and above Ep. The first one 
is mostly P p-derived while the other two are more hy- 
bridized with higher contributions from Co d-states. The 
energy alignments of these bands at T point are different 
for different Co derived skutterudites passing from (1)- 
(3)-(2) in CoSb 3 to (3)-(2)-(l) in CoAs 3 to (2)-(l)-(3) in 
CoP 3 . 

What we need to "realize" a favorable band structure 
for p-type TE properties from C0P3 is to lower the sin- 
gle degenerate band (1) until it crosses band (2) and 
arrives near to the heavy mass bands forming the bot- 
tom of the pseudogap so that we obtain a semiconductor 
with potentially high Seebeck coefficients analogous to 
La(Fe,Co)^Sbi2. The effect of gap opening by means 
of filling that provides TE performance for antimonides 
relies on the interaction between rare earth /-states or res- 
onances and the crossing band. This effect is strong be- 
cause the wavefunctions of the pseudogap crossing band 
have /-like symmetry as discussed in detail in Ref. [l3| 

Our results on LaF e^Pn show that the crossing band 
is indeed pushed down, as expected, by repulsion of La 
/-resonance states at about 3.0 eV above Ep (Fig J2| lower 
panel) resulting in a small direct gap between bands 
(1) and (2): E g — 98 meV. However the top of the 
band crossed by the Fermi level is not close to any 
lower heavy bands (Fig||) so that p-type thermopower 
can not be high enough at a reasonable band filling for 
La(Fe, C6)^Pyi assuming rigid band behavior upon al- 
loying with Co as found in La(Fe, Co)iSbi2 (N.B. strong 
non-rigid band behavior, which we do not expect, would 
also be detrimental to TE as it would indicate strong al- 
loy scattering and low carrier mobility). A study of the 
effective masses (Tab.|) for the three bands closest to Ep 
points out also that because of the double degeneracy and 
the reasonably high m*, La{Fe,C6)^P\2 could be more 
interesting for n-type TE application, but only if ther- 
mal conductivity can be strongly reduced and well-filled 
n-type material with low defect concentrations and high 
mobility can be produced. 

The importance of four membered pnictogen rings 
for thermoelectricity and also for superconductivity^ 
suggests investigation of the A g Raman active phonon 
frequenciesEEl at Brillouin zone center that are associated 
to normal modes involving variations of the symmetry- 
independent parameters (u,v) in the skutterudites struc- 
ture or, in other words, distortion of the pnictogen rings. 
We obtain LDA structural parameters (uq = 0.1462 and 
vq = 0.3478) near the experimental ones and lo\ = 228 
cm -1 and LO2 = 177 cm -1 for C0P3. Similar calcula- 
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tions for LaFe^P^ give uq — 0.1537 and vq = 0.3522 
for LDA equilibrium parameters whereas uj\ = 189 cm -1 
and ui2 — 160 cm -1 . The difference presumably reflects 
La-P interactions. 

In summary we have presented electronic structure cal- 
culations for C0P3 and LaFe^Pn. These show that while 
C0P3 is a metal, a gap is opened up upon filling with La. 
Nonetheless the band structure does not allow for high 
p-type thermopowers with reasonable carrier concentra- 
tions. The conduction bands are more favorable having 
a degenerate heavy mass structure, though we note that 
n-type filled skutterudites are difficult to prepare. 

We are grateful for useful discussions with R. S. Feigel- 
son. This work is supported by ONR and DARPA. 
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TABLE I. Effective masses for LaFaPw along the high 
symmetry directions in BZ. E(F) denotes the band energy (in 
eV) at T relative to E F (see Fig]|). 



FIG. 1. Band structure for C0P3 in an energy window cen- 
tered at Ef. Notice the metallic character due to band cross- 
ing at r. 
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FIG. 2. DOS for C0P3 (upper panel) and for LaFe 4 P 12 
(lower panel), the energy is referred to the respective Fermi 
level. Fhe total DOS (solid line) and the projections on 
different relevant atomic components that are sketched: P 
p-components (light gray shadow), Co d-components (dashed 
line) and resonant La /-components (dark gray shadow) . 

FIG. 3. Band structure for LaFdPw in an energy window 
including Ef- Notice the small gap opened at V. 
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